Nesting properties and anisotropy of the Fermi surface of LuNi 2 B 2 C 
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The rare earth nickel borocarbides, with the generic for- 
mula _RNi2B2C, have recently been shown to display a rich 
variety of phenomena. Most striking has been the competition 
between, and even coexistence of, antiferromagnetism and su- 
perconductivity. We have measured the Fermi surface (FS) 
of LuNi2B2C, and shown that it possesses nesting features 
capable of explaining some of the phenomena experimentally 
observed. In particular, it had previously been conjectured 
that a particular sheet of FS is responsible for the modulated 
magnetic structures manifest in some of the series. We report 
the first direct experimental observation of this sheet. 



The discovery of superconductivity in the rare earth 
nickel borocarbides [[l] has raised fundamental questions 
regarding its nature. The presence of moment-bearing 
rare-earth atoms in addition to high Ni contents makes 
the very existence of superconductivity surprising, and 
some of the phenomena associated with its interplay 
with magnetism have not been observed in any other 
superconducting material. With the generic formula 
.RN12B2C, some members of the system (R —Y, Lu, Tm, 
Er, Ho and Dy), exhibit moderately high values of T c 
(~16K, for R— Y,Lu), and are either antiferromagnetic 
or non-magnetic at low temperature. Two other systems, 
with R =Gd, Tb, are not superconducting, while the Yb 
compound shows heavy fermion behavior 

One of the most striking observations has been the an- 
tiferromagnetic ordering, and its competition (and even 
coexistence) with superconductivity, conjectured to be 
driven by a nesting feature in the Fermi surface (FS) 
||. Neutron and x-ray scattering techniques have re- 
vealed incommensurate magnetic structures in supercon- 
ducting Er |Q and Ho || compounds, and in the non- 
superconducting Tb H and Gd compounds, charac- 
terized by a wave vector Q m ~ (0.55,0,0). In the Lu 
compound, the 4/ band is fully occupied and therefore 
it is non-magnetic. Since the / electrons occupy local- 
ized core-like states, the FS is expected to be similar 
to that of the other compounds 0; unfettered by the 



complications introduced by the presence of magnetism, 
the Lu compound is ideal for investigating the origin of 
the superconductivity in these materials. Where there 
is magnetic order, the mutual influence of the moments 
must occur through indirect Ruderman-Kittel-Kasuya- 
Yosida (RKKY) type interactions; the resulting mag- 
netic structures would then be determined by maxima 
in the generalized magnetic susceptibility of the conduc- 
tion electrons. Band theoretical calculations of this sus- 
ceptibility, x(q)j m LUM2B2C show a peak around the 
aforementioned vector, Qm, indicating that the magnetic 
ordering in those other compounds is a result of a com- 
mon FS nesting feature ||. Moreover, the presence of 
strong Kohn anomalies in the phonon dispersion curves 
of LuNi2B2C for wave vectors close to Qm has lent ad- 
ditional support to this conjecture ||. Interest in the 
FS topology has been enhanced by the observation of a 
four- fold symmetry in the anisotropy of the upper critical 
field of LUM2B2C |j] . This has been interpreted in terms 
of (a) the anisotropy of the underlying FS topology, and 
(b) the presence of possible three-dimensional d-wave su- 
perconductivity JTo| . In the absence of any corroborative 
evidence for the latter suggestion, we shall concentrate 
on the former. The appearance of the square flux-line 
lattice (FLL) has been successfully explained through 
nonlocal corrections to the London theory [^TJ. In gen- 
eral, there is a nonlocal relationship between the super- 
current, j, and the vector potential, A, of the magnetic 
field in a superconductor, arising from the spatial ex- 
tent (~ £0) °f the Cooper pair [fl2|| . In this scenario, 
the shape of the FS (which influences the coupling of the 
supercurrents with the crystal lattice) is responsible for 
the square FLL observed at high fields in the mixed state 
and which reverts to the triangular FLL at lower fields 
p3| . However, despite extensive theoretical work based 
on the band structures of the non-magnetic Y and Lu 
compounds 14-16.3], very few experimental determina- 
tions of the electronic structure exist, and many of the 
predictions, including a "nestable" FS sheet, remain un- 
verified. Most importantly, the de Haas-van Alphen ex- 
periments performed on YNI2B2C in the superconducting 
p7[ and normal states |18| have not delivered the topo- 
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logical information necessary to isolate any such features 
in the FS. 

In this Letter, we present a joint experimental and the- 
oretical study of the FS of LuNi2B2C. The experiment re- 
veals the first direct evidence for the presence of a sheet 
capable of nesting f|. Furthermore, we will show that 
the anisotropy of this FS, as determined by our band 
structure calculations, is consistent with Kogan's model 
for the FLL (ll). 

The occupied momentum states, and hence the FS, 
can be accessed via the momentum distribution using the 
2-Dimensional Angular Correlation of electron-positron 
Annihilation Radiation (2D-ACAR) technique A 
2D-ACAR measurement yields a 2D projection (integra- 
tion over one dimension) of an underlying two-photon 
momentum density, 



p(p) = « J (k)|C Gj (k)| 2 <5(p-k-G) 



(1) 



j.k.G 



where v? (k) is the electron occupation density in k-space 
in the j^ 1 band, the CG,j(k) are the Fourier coefficients 
of the interacting electron-positron wave function prod- 
uct and the delta function expresses the conservation of 
crystal momentum. p(p) is a single-centered distribution 
having the full point symmetry of the the crystal lattice 
in question. In a metal, this distribution contains dis- 
continuities at various points pp = (kp + G) when a 
band crosses the Fermi level, £p. When the FS is of 
paramount interest, the Lock- Crisp- West procedure |2C| l 
is often followed. Here the various FS discontinuities are 
superimposed by folding p(p) (or its measured projec- 
tions) back into the first Brillouin zone (BZ). The re- 
sult is a new k-space density, J2j k ^ s ) Sg |CG,j(k)| 2 , 
which aside from the factor ^ G |CG.j(k)| 2 (usually a 
weak function of k within each band) is simply the elec- 
tron occupation density. This well-established technique 
has recently been used by some of the present authors 
to identify and measure the FS topology in pure Y 
and in disordered Gd-Y alloys |p2[| . The virtue of the 
2D-ACAR technique in such studies is that it reveals di- 
rectly the shape of the FS, and hence any propensity for 
nesting. 

The experiments were performed on a single crystal 
of LuNi 2 B 2 C, grown by a Ni 2 B flux method Q. The 
sample was cooled to ^50K, at which temperature the 
overall momentum resolution of the Bristol 2D-ACAR 
spectrometer corresponded to ~10% of the larger BZ di- 
mension. Projections were measured along two different 
crystallographic directions, a and c. More than 400 mil- 
lion (effective) counts were collected in each spectrum. 

The spin-dependent momentum densities were calcu- 
lated using the linearized muffin-tin orbital (LMTO) 
method within the atomic sphere approximation (ASA), 
including combined-correction terms p5| . The exchange- 
correlation part of the potential was described in the 



local density approximation. The self-consistent band- 
structure was calculated at 594 fc-points in the irreducible 
l/lQ th part of the BZ, which to simplify the calculations 
had the same volume as the standard BZ but a simpler 
tetragonal shape. We used a basis set of s, p, d and / 
functions for the Lu, and a reduced (s,p,d) basis set for 
Ni, B and C. For the calculation of the density-of-states 
(DOS), and the FS, a denser mesh of 3696 fc-points was 
used. The lattice parameters used were experimental val- 
ues of a = 3.464 A and c = 10.631 A |||. Electronic 
wave functions from the dense fc-mesh were then used to 
generate the electron-positron momentum densities, us- 
ing 1149 reciprocal lattice vectors. A full description of 
the technique is given in Refs. J2^| and 27|. 

The band-structure obtained (not shown) is very simi- 
lar to those of Pickett and Singh jl6) and Mattheiss jlJJ , 
and shows that the electronic structure is certainly three- 
dimensional (compared to the two-dimensional features 
exhibited in the high-T c cuprates, e.g. |p8| ); experimen- 
tally, one observes an almost isotropic resistivity [p9| , 
which is consistent with this picture. These bands pre- 
dict a rather complicated FS (shown in Fig. [l]), the prin- 
cipal character being Lu d, with some Ni d, and B and 
C p character The first sheet is a very small 

electron pocket centered at T, while the second is slighly 
larger and rather more cuboid in shape. It is the third 
sheet that possesses the nesting properties previously re- 
marked upon H , but it can also be seen that the second 
sheet has a square cross section between M and R. 

In Figs. H and || the BZ electron occupancies (both 
calculated and experimental) are shown, projected along 
the c and a axes, respectively. The occupied states are 
indicated by the white areas, and unoccupied by black. 
The agreement between the positron experiment and the 
LMTO theory is excellent; one can clearly discern the 
presence of the T-centered electron sheets. The most 
striking features of Fig. || are the electron surfaces with 
square cross-sections at the corners of the projected BZ. 
Also noteworthy are the shapes described by the con- 
tour lines around this sheet; the protuberances point 
towards the T-point in both the calculation and exper- 
iment. Fig. H shows fewer features, but demonstrates 
that the agreement persists when the projection is along 
a different direction. 

Considerable investment has been made in establish- 
ing methods of reliably and accurately locating the FS 
in 2D-ACAR data (30). Recently, it has been demon- 
strated that it is possible to accurately caliper the FS 
from 2D-ACAR data using a criterion based upon edge- 
detection algorithms employing band-pass filtering or 
Maximum Entropy techniques pl]-|23]|. In Fi g. fl| we 
present the experimental FS thus derived [ pl| , |22| ], to- 
gether with a (001) section (passing through the r-point) 
of the calculated third FS sheet. The nesting is indi- 
cated by the arrow, and is calipered from the experi- 
mental data to be 0.54 ± 0.02 x {2-k/o) (our calculation 
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gives 0.56 x (2ir/a)). Given that the experimental data 
represent a projection of the FS, one does not expect 
a perfect match between the top and bottom halves of 
the figure. However, the nesting feature is unequivo- 
cally revealed, and its size is in excellent agreement with 
the pronounced Kohn anomalies close to a wave vector 
Qni ~ (0.5,0,0) ||], and the observed incommensurate 
ordering with Q m « (0.55,0,0) found in the antiferro- 
magnetic compounds. In conjunction with the calcula- 
tion, it was possible to estimate from the experiment 
that the fraction of the FS that would be able to par- 
ticipate in nesting is (4.4 ± 0.5) %, which is consistent 
with the small increase in the resistivity observed when 
the current is along [100] 

Having shown that the calculated FS topology is in 
good agreement with our experiment, we derived perti- 
nent quantities from our band structure. The DOS at the 
Fermi energy, N(Ep) — 4.3 (eV cell) -1 , is similar to val- 
ues presented by Pickett and Singh |l(J and Mattheiss 
0] (4.8 (eV cell)" 1 ). Most importantly, some of the 
quantities essential to Kogan's theory |ll[] can be derived 
from the band structure itself. For the calculation of 
the Fermi velocities, a special mesh was used, compris- 
ing six additional points around each original fc-point, 
to enable accurate evaluation of the relevant derivatives. 
The Fermi velocities in the ab plane and in the c direc- 
tion [|§ were {v 2 Fab ) 1/2 = 2.62 x 10 7 cm s" 1 ,^) 1 / 2 = 
2.23 x 10 7 cm s , with masses m a (, = 0.91, m c = 1.26, 
implying an average out-of-plane anisotropy in the upper 
critical field ((H< 100> +H< 110> )/2H< 001> ) of 1.17 com- 
pared to an experimental value of 1.16 A calculation 
for the isovalent compound YNi 2 B 2 C predicted 1.02 for 
this anisotropy, in good agreement with the experimental 
value of 1.01 |33| and highlighting its sensitivity to the 
band structure. 

In conclusion, we have shown experimentally that the 
FS topology of LuNi 2 B 2 C does support nesting, thereby 
accounting for the anomalies observed in its phonon spec- 
trum ||] , and the propensity for magnetic ordering found 
in some of the other rare-earth nickel borocarbides Q . In 
addition, our own calculations of the electronic structure 
yield a FS in excellent agreement with the experiments, 
whose anisotropy is consistent with the observation of a 
square FLL, and with the observed behavior of the upper 
critical field ||. 
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FIG. 2. Experimental (left) and calculated electron density 
(right) projected along the [001] direction, with the symmetry 
points of the simple tetragonal BZ indicating the projection 
path. Black signifies holes, and white represents electrons. 
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FIG. 3. Experimental (left) and calculated electron density 
(right) projected along the [010] direction. 
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FIG. 1. The three sheets of the FS of LuNi 2 B 2 C, shown 
within the simple tetragonal BZ (see text), with its usual 
symmetry points labelled. 
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FIG. 4. The experimental (top) and calculated (bottom) 
FS topology of LuNi2B2C. The calculation is of the FS in the 
third band in the (001) plane through the r-point. The arrow 
indicates the nesting feature. 
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